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Abstract. We study the reliability of dark-matter halo detections with three different linear filters applied to weak-lensing data. 
We use ray-tracing in the multiple lens-plane approximation through a large cosmological simulation to construct realizations 
of cosmic lensing by large-scale structures between redshifts zero and two. We apply the filters mentioned above to detect 
peaks in the weak-lensing signal and compare them with the true population of dark matter halos present in the simulation. We 
confirm the stability and performance of a filter optimised for suppressing the contamination by large-scale structure. It allows 
the reliable detection of dark-matter halos with masses above a few times 10'^ Mq with a fraction of spurious detections 
below ~ 10%. For sources at redshift two, 50% of the halos more massive than ~ 7x 10" h^' Mq are detected, and completeness 
is reached at ~ 2 x lO'"* Mr... 



1. Introduction 

How reliably can dark-matter halos be detected by means of 
weak lensing, and what selection function in terms of mass and 
redshift can be expected? This question is important in the con- 
text of the analysis of current and upcoming wide-field weak- 
lensing surveys. This subject touches upon a number of scien- 
tific questions, in particular as to how the non-linear growth 
of sufficiently massive structures proceeds throughout cosmic 
history, whether galaxy-cluster detection based on gas physics 
agrees with or differs from lensing-based detection, whether 
dark-matter concentrations exist which emit substantially less 
light than usual or none at all, what cosmological information 
can be obtained by counting dark-matter halos, and so forth. 

As surveys proceed or approach which cover substantial 
fractions of the sky, such as the CFHTLS survey, the upcom- 
ing Pan-STARilS surveys, or the planned surveys with the 
DUNE or SNAP satellites, automatic searches for dark-matter 
halos w ill r outinely be c a rried o ut, see for example lErben et all 
(l2000h and lErben et al.l (l2003h It is important to study what 
they are expected to find. 

Several different methods for identifying dark-matter ha- 
los in weak-lensing data have been proposed in recent years. 
They can all be considered as variants of linear filtering tech- 
niques with different kernel functions. Particular examples are 
the aperture mass with th e radial filter functions proposed 



by Schneider et al. ( 1998b and modifi ed by Schirmer et al.. 



(|2004 and iHennawi & Spergell (l2005b . and the filter opti- 



embedded in (IMaturi et al.L l2005h . The non-negligible con- 
tamination by the large-sca le structure was alread y note d 
by RebUnskv & BartelmannI tl999) and IWhite et all (l2002h . 



and Hoekstra 



(1200 II) quantified its impact on weak-lensing 
mass determinations. iHennawi & Spergell (l2005b showed that 
the redshift of background galaxies can be used to improve 
the number of reliable detections. An approach alternative to 
match ed filters is based on the peak statistics of convergence 



maps (IJain & Van Waerbeke , l2000l). e.g. obtained with the 



Kaiser-Squires inversion techiiique (IKaiser & Squires L 11993 



Kaiser. Squires & Broadhurst L 119951) or variants thereof. 



In this paper, we evaluate three halo-detection filters in 
terms of their performance on simulated large-scale structure 
data in which the dark-matter halos are of course known. One 
of the filters is specifi cally designed to o ptimally suppress 



the LSS contamination ( Maturi et al.L 12005 ). This allows us to 



mised for separating the weak-lensing signal of dark-matter 
halos from that of the large-scale structures (LSS) they are 



quantify the completeness of the resulting halo catalogues, the 
fraction of spurious detections they contain, and the halo se- 
lection function they achieve. In particular, we compare the 
performance of the three filters mentioned in order to test and 
compare their refiabifity under a variety of conditions. 

We summarise the required aspects of lensing theory in 
Sect. 2 and describe the numerical simulation in Sect. 3. The 
weak-lensing filters are discussed in Sect. 4, and results are 
presented in Sect. 5. We compare suitably adapted simulation 
results to the GaBoDS data in Sect. 6, and present our conclu- 
sions in Sect. 7. 
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2. Lensing theory 



2.2. Multiple lens-plane theory 



In this section, we briefly summarise those aspects of gravita- 
tional lensing that are relevant for the present study. For a more 
detailed discussion on the the ory of gravitational lensing we re - 
fer the reader to the review bv lBartelmann & Schneiden (1200 lb . 



2.1. Lensing by a single plane of matter 

We start with the deflection of light rays by thin structures 
in the universe. This thin-screen approximation applies when 
the physical size of the gravitational lens is small compared to 
the distances between the observer and the lens and between 
the lens and the sources. Accordingly, we project the three- 
dimensional matter distribution of the lens p{0, w) on the lens 
plane and obtain the surface mass density 



2(0) 



p(6, w)dw. 



(1) 



where is the angular position vector on the lens plane and w 
is the radial coordinate distance along the line-of-sight. 
The critical surface density is defined as 



4nG DdD, 



(2) 



where Ds, Dd and Dds are the angular-diameter distances be- 
tween the observer and the source, between the observer and 
the lens, and between the lens and the source, respectively. 
Rescaling the surface density with the critical surface density 
we obtain the convergence 



(3) 



A thin lens is fully described by the lensing potential 
which is related to the convergence through the two- 
dimensional Poisson equation 



2k. 



(4) 



The deflection angle is the gradient of the lensing potential, 



(5) 



We now generalise the previous formalism to the case of a con- 
tinuous distribution of matter filling the volume between the 
observer and the sources. We can divide this volume into a 
sequence of equidistant sub-volumes whose depths along the 
line-of-sight are small compared to the cosmological distances 
between the observer and the centres of the sub-volumes, and 
between those and the sources. Doing so, the thin screen ap- 
proximation remains locally valid. If the light cone is narrow 
enough, the matter distribution within each sub-volume can be 
replaced by a two-dimensional matter sheet. Thus, light rays 
are approximated by a sequence of straight lines between the 
planes between the observer and the sources, where is the 
number of sub-volumes into which the space between observer 
and sources has been divided. In the following discussion, we 
assume that the sources lie on the (A^ + l)-th plane. 

The deflection angles on each plane can be computed 
by spatiall y differencing the corres ponding lensing potential. 
Following IHamana & Mellierl (120011) . the gravitational poten- 
tial of the matter within each sub-volume is decomposed into 
a background and a perturbing potential. The equation relating 
the lensing potential to the mass distribution responsible for 
lensing on the i-th plane is very similar to Eq.|4] but the surface 
density is substituted by the density contrast of the projected 
matter 



srProj 



(8) 



where p, is the three-dimensional density field in the i-th sub- 
volume, p is the mean density of the universe, and 6, - {pi - 
p)lp is the density contrast of the three-dimensional mass dis- 
tribution in the i-th sub- volume. The Poisson equation (HJi is 
thus 



VV/(0) = 



%nGp ^ 



3c2 



\0) = 30 



\0) 



(9) 



and the deflection angle is 
or/ = V(A/ . 



(10) 



A light ray propagating from the source to the observer is 
thus deflected on each plane by the amount aiiOi), where 0, 
indicates the position where the light ray intercepts the /-th lens 
plane. The lens equation, relating the positions of a light ray on 
the A^-th and the first planes can be constructed iteratively by 
summing the deflections on all intermediate planes. It becomes 



The point on the lens plane is mapped onto the point P on 
the source plane given by the lens equation. 



P^e- a(0). 



(6) 



where a(ff) - Ddsl Dsa{0) is the reduced deflection angle. 

The complex shear, y = (yi -H iyi), is also obtained from 
the lensing potential. Its components ji and 72 are 



1 



del J 



72 



(7) 



. 1 J-^n'-^is 



(11) 



which represents the generalisation of Eq.|6]to the multi -plane 
case. In the previous equation, we introduced D,„ and D,s to 
indicate the angular diameter distances between the i-th and 
the A^-th planes, and between the i-th plane and the sources, 
respectively. 

The lens mapping between the i-th and the first planes is 
described by the Jacobian matrix 

d0i 



Ai = 



30, 



(12) 
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The matrix describing the mapping through a sequence of 
planes is obtained by recursion. It is given by 



An - I - > UiA 



where 



(13) 



(14) 



On the source plane, we define the Jacobian matrix A^r+i by 



1-K- 

-72 - 



■ 7i -72 + w 
Li) 1 - /<■ + 7l 



(15) 



This is not necessarily symmetric since it is the product of two 
symmetric matrices. The asymmetry is given by the rotation 
term lj. This term only appears in the multiple lens-plane the- 
ory because the Jacobian matrix is symmetric for a single lens 
plane. The terms k and y appearing in Eq.[T5]are now the effec- 
tive convergence and the effective shear, respectively. 

3. The numerical simulation 

3.1. The cosmological box 

The cosmological simulation used in this study is the result of a 
hydrodynami cal, A^-body si mulation, carried out with the code 
GADGET-2 (Springe!, 20Q5V It has been described and used in 



several previous studies ( Murante et al.Ll2004l:lRoncarelli et al. 



20061) . We only briefly summarise here some of its characteris- 
tics. A more detailed discussion can be found in the paper by 
Borganiet al.i ( .2004) . 

The simulation represents a concordance ACDM model, 
with matter density parameter f2,„ = 0.3 and a contribution 
from the cosmological constant Qa = 0.7. The Hubble parame- 
ter is h = i/o/100 - 0.7 and a baryon density parameter Qbar = 
0.04 is assumed. The normalisation of the power spectrum of 
the initial density fluctuations, given in terms of the rms density 
fluctuations in spheres of 8 /z~'Mpc, is erg - 0.8, in agreement 
with the most recent constraints from weak lensing and from 
the observations of the Cosmic Microwa ve Background (e.g. 



Hoekstra et all 120061; iSpergel et al.L 120061) 



The simulated box is a cube with a side length of 
192/i"'Mpc. It contains 480^ particles of dark matter and an 
equivalent number of gas particles. The Plummer-equivalent 
gravitational softening is set to epi - 7.5kpc/h comoving be- 
tween redshifts two and zero, and chosen fixed in physical units 
at higher redshift. 

The evolution of the gas component is studied including 
radiative cooling, star formation and supernova feedback, as- 
suming zero metalicity. The treatment of radiative cooling as- 
sumes an optically thin gas composed of 76% hydrogen and 
24% of helium by mass, plus a time-dependent, photoionis- 
ing uniform UV background given by quasars reionising the 
Universe at z a; 6. Star formation is implemented using the hy- 
bri d multiphase model for the i nterstellar medium introduced 
by ISpringel & HernquistI (120031) . according to which the ISM 
is parameterised as a two-phase fluid consisting of cold clouds 
and hot medium. 



The mass resolution is 6.6 x IO'^Mq/Zi for the cold dark 
matter particles, and 8.9 x IQ^Mq/H for the gas particles. This 
allows resolving halos of mass 10'^/;"'Mq with several thou- 
sands of particles. 

Several snapshots are obtained from the simulation at scale 
factors which are logarithmically equidistant between fl,„,- - 
0.1 and afi„ - 1. Such snapshots are used to construct light- 
cones for the following ray-tracing analysis. 

3.2. Construction of the light-cones 

Aiming at studying light propagation through an inhomoge- 
neous universe, we construct light-cones by stacking snapshots 
of our cosmological simulation at different redshifts. Each 
snapshot consists of a cubic volume containing one realization 
of the matter distribution in the ACDM model at a given red- 
shift. However, since they are all obtained from the same initial 
conditions, these volumes contain the same cosmic structures 
in diff'erent stages of their evolution. Such structures are ap- 
proximately at the same positions in each box. Hence, if we 
want to stack snapshots in order to build a light-cone encom- 
passing the matter distribution of the universe between an ini- 
tial and a final redshift, we cannot simply create a sequence of 
consecutive snapshots. Instead, they must be randomly rotated 
and shifted in order to avoid repetitions of the same cosmic 
structures along one line-of-sight. This is achieved by apply- 
ing transformations to the coordinates of the particles in each 
cube. When doing so, we consider periodic boundary condi- 
tions such that a particle exiting the cube on one side re-enters 
on the opposite side. 

One additional problem in stacking the cubes is caused by 
the fact that, as they were written at logarithmically spaced 
scale factors, consecutive snapshots overlap with each other 
by up to two-thirds of their comoving side-length (at the low- 
est redshift). Thus, we have to make sure to count the matter 
in the overlapping regions only once. For doing so, we chose 
to remove particles from the later snapshot. The choice of the 
particles to remove from the light cone is not critical, since 
snapshots are relatively close in cosmic time. Several tests have 
confirmed this expectation. 

Hence, the light-cone to a given source redshift Zs is con- 
structed by filling the space between the observer and the 
sources with a sequence of randomly rotated and shifted vol- 
umes. As explained in Sect. 12.21 if the size of the volumes is 
small enough, we can approximate the three-dimensional mass 
distribution in each volume by a two-dimensional mass distri- 
bution. This is done by projecting the particle positions on the 
mid-plane through each volume perpendicular to the line-of- 
sight. Such planes will be used as lens planes in the following 
ray-tracing simulations. 

The opening angle of the light-cone is defined by the an- 
gle subtending the physical side-length of the last plane before 
the source plane. For sources at z, = 1 and z, = 2, this corre- 
sponds to opening angles of 4.9 and 3.1 degrees, respectively. 
In principle, tiling snapshots at constant cosmic time allows the 
creation light-cones of arbitrary opening angles. However, this 
is not necessary for the purposes of the present study. 
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Fig. 1. Sketch illustrating the construction of the light cones. 
A sequence of lens planes (vertical lines) is used to fill the 
space between the observer (O) and the sources on the (A^ + 1)- 
th plane. The aperture of the light cone depends on the distance 
to the last lens plane. At low redshifts, only a small fraction of 
the lens planes enters the light-cone (dark-gray shaded region). 
This fraction increases by reducing the redshift of the sources, 
increasing the aperture of the light cone (light-gray shaded re- 
gion). 



If the size of the light-cone is given by the last lens plane, 
increasingly smaller fractions of the remaining lens planes will 
enter the light cone as it approaches z = (a = 1, see FigdJ. 

3.3. Halo catalogues 

Each simulation box contains a large number of dark-matter 
halos. For our analysis, it is fundamental to know the location 
of the halos as well as some of their properties, such as their 
masses and virial radii. Thus, we construct a catalog of halos 
for each snapshot. The procedure is as follows. We first run a 
friends-of-friends algorithm to identify the particles belonging 
to a same group. The chosen linking length is 0.15 times the 
mean particle separation. Then, within each group of linked 
particles, we identify the particle with the smallest value of the 
gravitational potential. This is taken to be the centre of the halo. 
Finally, we calculate the matter overdensity in spheres around 
the halo centre and measure the radius that enclosing an aver- 
age density equal to the virial density for the adopted cosmo- 
logical model, pvir = ^c(z)pcm(z), where Pcrit(z) is the critical 
density of the universe at re dshift z, and the overdensity Adz) 
is calculated as described in Eke et al.l (l200lb . 

We end up with a catalogue containing the positions, the 
virial masses and radii, and the redshifts of all halos in each 
snapshot. The positions are given in comoving units in the co- 
ordinate system of the numerical simulation. They are rotated 
and shifted in the same way as the particles during the construc- 
tion of the light-cone. The positions of the halos in the cone are 
finally projected on the corresponding lens plane. 

In Fig|2l we show the mass functions of the dark-matter ha- 
los normalised to one square degree and contained in the light- 
cones corresponding to Zs = 1 (solid line) and Zs = 2 (dashed 
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Fig. 2. Number of halos per mass bin per square degree. The red 
and green curves show the halo mass distribution for sources at 
Zs = 1 and Zs = 2, respectively. 

line). Obviously the light-cones contain a large number of low- 
mass halos (~ 10^' - IO'^/i^'Mq) which are expected to be 
undetectable through weak lensing. On the other hand, a much 
lower number of halos with mass M > are po- 

tential lenses. We note that the numbers of haloes with masses 
larger than 5 x lO'^M© are approximately equal in both light 
cones, because such haloes are mainly contained in the low- 
redshift portion of the volume which is common to both light 
cones. 

We ignore the intracluster gas here because it contributes 
about one order of magnitude less mass than the dark mat- 
ter and therefore does not significantly affect the weak-lensing 
quantities. 

3.4. Ray-tracing simulations 

The lensing simulations are carried out using standard ray- 
tracing techniques. Briefly, starting from the observer, we trace 
a bundle of 2048 x 2048 light rays through a regular grid cover- 
ing the first lens plane. Then, we follow the light paths towards 
the sources, taking the deflections on each lens plane into ac- 
count. 

In order to calculate the deflection angles, we proceed as 
follows. On each lens plane, the particle positions are interpo- 
lated on regular grids of 2048 x 2048 cells using the tri a ngular - 



shaped-cloud (TSC) scheme dHocknev & Eastwood 119881) . 
This allows to avoid sudden discontinuities in the lensing 
mass distributions, that would lead to anomalous deflections 



of the light rays dMeneghetti et al. ■ I2OOOI: iHamana & Mellieit 



2OOII) . The resulting projected mass maps, M', , where /, m 



1, 2048 and / = 1, A^, are then converted into maps of the 
projected density contrast. 



M' 

6'!""'' = —-Li 
'"' Aip 



(16) 



where A, and L, are the area of the grid cells on the i-th plane 
and the depth of the i-th volume used to build the light cone, 
respectively. 
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The lensing potential at each grid point, i/'J^^^, is then calcu- 
lated using Eq. |9l Owing to the periodic boundary conditions 
of the density-contrast maps, this is easily solvable using fast- 
Fourier techniques. Indeed, Eq. |9] becomes linear in Fourier 
space. 



(17) 



where k is the wave vector and and 6'"'°^ are the Fourier 
transforms of the lensing potential and of the projected density 
contrast, respectively. Using finite-differencing schemes, we fi- 
nally obtain rnaps of the deflection angles on each plane, a'^^^^ 



dPremadi et 31.11199 81) 



The arrival position of each light ray on the source plane is 
computed using Eq. [TT] which incorporates the deflections on 
all preceding lens planes. However, the ray path intercepts 
the lens plane at arbitrary points, while the deflection angles 
are known on regular grids. Thus, the deflection angles at the 
ray position are calculated by bi-linear interpolation of the de- 
flection angle maps. 

Again using finite differencing schemes, we employ 
Eqs.[T2]to[T5]to obtain maps of the eflTective convergence and 
shear. 



3.5. Testing the ray-tracing code 

We test the reliability of the ray-tracing code by comparing the 
statistical properties of several ray-tracing simulations with the 
theoretical expectations for a ACDM cosmology. In these tests, 
we assume that all source redshifts are Zs = 1.5. For this source 
redshift, the light cone spans a solid angle of roughly 3.6^ - 13 
square degrees on the sky. We perform ray-tracing through 60 
diflerent light-cones in total. 

In Fig. [3] we show the power spectra of the eff'ective conver- 
gence and the shear, obtained by averaging over all different re- 
alizations of the light-cone. These are given by the solid and by 
the dotted lines, respectively. The theoretically expected power 
spectrum is shown as the dashed line. As expected, the conver- 
gence and the shear power spectra are equal. We note that they 
agree with the theoretical expectation over a limited range of 
wave numbers. Indeed, they deviate from the theoretical power 
spectrum for k < 200 and for k > 20000. These two values 
of the wave vector define the reliability range of these simula- 
tions and are both determined by numerical issues. On angular 
scales > 1°, we miss power because of the small size of the 
simulation box, while on angular scales smaller than < 1' we 
suff'er from resolution problems due to the finite resolution of 
the ray and the mass grids. 



3.6. Lensing of distant galaxies 

Using the effective convergence and shear maps obtained from 
the ray-tracing simulations, we are now able to apply the lens- 
ing distortion to the images of a population of background 
sources. 
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Fig. 3. Numerical power spectra of the effective convergence 
(solid line) and of the shear (dotted line) obtained by averag- 
ing over 60 different light cones corresponding to a solid angle 
of ~ 13 square degrees. The power spectrum expected for a 
ACDM model with the same cosmological parameters as the 
simulation is given by the dashed line. The errorbars are shown 
only for the effective convergence power spectrum, but are of 
equivalent size for the shear power spectrum. 

In the weak-lensing regime, a sufficiently small source with 
intrinsic ellipticity ej is imaged to have an ellipticity 



es +; 



(18) 



where g is the complex reduced shear, and the asterisk denotes 
complex conjugation. We adopt here the standard definition of 
ellipticity, e - ei + iei - {a - b)/(a + b)e'^'f', were a and b are 
the semi-major and the semi-minor axes of an ellipse fitting the 
object's surface-brightness distribution. The position angle of 
the ellipse's major axis is (p. 

The complex reduced shear is defined as 



7(0) 



(19) 



In the weak-lensing regime, /c ^ 1 and y ~ g- Equation [18] 
illustrates that the lensing distortion is determined by the re- 
duced shear. Since the background galaxies are expected to be 
randomly oriented, the expectation value of the intrinsic source 
ellipticity ej is assumed to vanish. Thus, the ellipticity induced 
by lensing can be measured by averaging over a sufficiently 
large number of galaxies and the expectation value of the ellip- 
ticity is equal to the reduced shear, {e) - g. 

In order to generate a mock catalogue of lensed sources, 
galaxies are randomly placed and oriented on the source plane. 
Their intrinsic ellipticities are drawn from the distribution 



P(\es\) 



expKl 



;r(r2|exp(l/cr2)-l| 



(20) 



where cr^^ = 0.3. We assume a background galaxy number 
density of Ug - 30 arcmin 2. Observed elUpticities are ob- 
tained from Eq. [T8]by interpolating the effective convergence 
and shear at the galaxy positions. This procedure results in 
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catalogues of lensed galaxies for each source redshift chosen, 
where galaxy positions and ellipticities are stored. 

4. Weak lensing estimators 

We investigate the performances of three weak-lensing estima- 
tors which have been used so far for detecting dark-matter ha- 
los through weak lensing. These are the c lassical aperture mass 
dSch neider. 1996; Schneide r et al. . 1998 ). an optimised version 
of it (.Schirmer et al.. .2004i) . and the receritly de veloped, opti- 
mal weak-lensing halo filter (IMaturi et al.ll2005h . More details 
on these three estimators are given below. 

All of them measure the amplitude of the lensing sig- 
nal A within circular apertures of size around a centre 0. 
Generalisations are possible to apertures of different shapes. In 
general, A is expressed by a weighted integral of the tangential 
component of the shear relative to the point 6, y,. The weight 
is provided by a filter function *P, such that 



A(6I,0)= A^ffyAe'^efViW -0\) 



(21) 



and the integral extends over the chosen aperture. The variance 
of the weak-lensing estimator is given by 



cr 



(2;r)2 / 



(22) 



where *l'(^) is the Fourier transform of the filter, and Pnik) the 
power spectrum of the noise. 

4.1. Aperture Mass 



The a perture mass was originally proposed by dSchneider , 
1996 V for measuring the projected mass of dark-matter con- 
centrations via weak lensing. It represents a weighted integral 
of the convergence. 



MAPTiO) = J d^e'K(0')u(\e' 



■o\) 



(23) 



The weight function U(d) is symmetric if the aperture is chosen 
to be circular, and it is compensated, i.e. 



Jo 



(24) 



Since the convergence is not an observable, the aperture 
mass is more conveniently written as a weighted integral of the 
tangential shear. 



Mapt{0) 



0)Tapt(|0' - e\) 



(25) 



where the function ^apt is related to the filter function U by 
the equation 



2 , , , 



) - U(0) . 



The variance cr], of Mapt is defined 



M, 



as 



^Mapt 



— d0'0''i' 

ng Jo 



■lpA\o'-e\) 



(26) 



(27) 



which takes into account the shot noise due to the fi- 
nit e number and the intrinsic ellipticities of the sources 



dBartelmann & Schneideiil200Ih . 

The shape of the filter function is usually chosen to have 
a compact support and to suppress the halo centre because the 
lensing measurements are more problematic there. Indeed, the 
weak-lensing approximation may break down and the cluster 
galaxies may prevent the ellipticity of background galaxies to 
be accurately measured. 



Schneider et alJ dl998l) propose the polynomial function 



(28) 



where H(x) is the Heaviside step function, and x - 0/0max 
is the radial angular coordinate in units of the radius, 0max, 
where ^'apt vanishes. Z is a free parameter which is usually 
set to Z = 1. Note that this filter function was designed es- 
pecially for measuring cosmic shear. However, se veral authors 



have used it for searches f or dark matter halos dErben et al 
20001: ISchirmer et alll2004 . 



More recently, other filter functions have been pro- 
posed which maximise the signal-to-noise ratio Mapt/ctmapt- 



Schneider et alJ (Il998h show that this is the c ase if Q mim 



ics the shear profile of the lens. For example, ISchirmer et al 



d2004f) propose a fitting formula that approxim ates the shear 



profil e of a Navarro-Frenk-White (NFW) halo (iNavarro et al 
19961) . Their filter function is 



tmhx/Xc 



X/Xc 



(29) 



whe re Xc is a parameter c ontrolling the shape of the filter (see 



also Padmanab han et al.L l2003t iHetterscheidt et all 120051) . In 



the rest of the paper, we will refer to this implementation of the 
ape rture mass as to the "optimi sed aperture mass". 



Hennawi & Spergel ( 2005^) included the photometric red- 



shifts of background sources, increasing the halo-detection sen- 
sitivity at higher redshifts and for smaller masses. Aiming at a 
comparison of different filters, we neglect this additional in- 
formation here. We can therefore not apply their tomographic 
approach, which is based on an NFW fitting formula. They also 
suggested using a Gaussian profile whi ch found application i n 



actual weak-lensing surveys (see e.g. Mivazaki et al. . 2002h 



but here we focus on the filter proposed by lMaturi et al. I d2005b 



whose shape is statistically and physically well motivated. 
4.2. Optimal Filter 

Maturi et al. I d2005h have recently proposed a weak-lensing fil- 
ter optimised for an unbiased detection of the tangential shear 
pattern of dark-matter halos. Unlike the optimised aperture 
mass, the shape of optimal filter is determined not only by the 
shear profile of the lens, but also by the properties of the noise 
affecting the weak lensing measurements. 

The measured data D is composed of the signal from the 
lens S and by the noise A^, and can be written as 



D{0) = 5(61) H- N(e) = At(0) + N(0) 



(30) 
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where A is the total amplitude of the tangential shear and 7(0) is 
its angular shape. The noise comprises several contributions 
that can be suitably modeled. 

The optimal filter accounts for the noise contributions be- 
cause it is constructed such as to satisfy two conditions. First, 
it has to be unbiased, i.e. the average error on the estimate of 
the lensing amplitude. 



0\) 



has to vanish: 



b = {A,,,-A)=A 



Second, the noise 



= 



a^ = {(A,„-Af)^b' + 



{2nf j 



mk)\^Pdk)d^k 



(31) 



(32) 



(33) 



has to be minimal with respect to the signal. 

The filter function *F satisfying these two conditions is 
found by combining them with a Lagrangian multiplier A. The 
variation L - cr^ + Ab is carried out, and the filter function ^' 
is found by minimising L. In Fourier space, the solution of this 
variational minimisation is 



'i'oPTik) = 



1 



\T(kW 
PN(k) ■ 



m_ 

PN(k) 



(34) 



where the hats denote the Fourier transform. The last equation 
shows that the shape of the optimal filter is determined by 
the shape of the signal, t, and by the power spectrum of the 



noise, P^- 



Maturi et al. 



(I2005h model the signal by assuming 
that clusters are on average axially symmetric and their 
shear profile resembles that of an NFW halo (see e.g . 
Bartelmannl 1 19961: IWright & Brainerdl llOOd iLi & Ostrikei . 
2002t iMeneghetti et all l2003h . Consequently, this filter is op- 
timised for searching for the same halo shape as the optimised 
aperture mass, even if the filter profile is different. 

The noise is assumed to be given by three contributions, 
namely the noise contributions from the finite number of back- 
ground sources, the noise from their intrinsic ellipticities and 
orientations, and the weak-lensing signal due to the large-scale 
structure of the universe. 

The first two sources of noise are characterised by the 
power spectrum 



PAk)^-^ 
1 n„ 



(35) 



which depends on the dispersion of the intrinsic ellipticities 
of the sources, cr^ , and on the number density of background 
galaxies, Ug. 

The statistical properties of the noise due to the lensing sig- 
nal from the large-scale structure of the universe are described 
by the power-spectrum of the eff'ective tangential shear This is 
related to the power-spectrum of the effective convergence by 



Py(.k) = ^PAk) . 



Thus, the total noise power spectrum is 
P^k) = PJk) + P,(k) , 



(37) 



where Py is determined by the linear theory of structure growth. 
Using the linear instead of the non-linear power spectrum 
avoids suppressing a substantial fraction of the signal from the 
non-linear structures we are searching for To further reduce 
any loss of signal in the filtering process, it would be possible 
to cut Py off at angular scales typical for galaxy clusters. Doing 
so, we found that this approach has a negligible impact on the 
final result. 



— 1 1 1 1 1 r 

OPT, excluding LSS 

OPT, suppressing LSS 

THS 

OAPT 

APT 

GAUSSIAN 




Fig. 4. Comparison of the different filter shapes used here and 
in the literature. The filter scales are those typically used in 
the literature. Note how the optimal filter (black solid curve) 
shrinks when the linear matter power spectrum is used to sup- 
press the LSS contributio n (red dashed curve). Interestingly, 
Hennawi & Spergell(l2005 l) found experimentally that the trun- 
cated NFW-shaped filter (cyan curve) performs best when 
scaled to the green curve (THS), which approximates the opti- 
mal filter (OPT, red curve) almost precisely. The advantage of 
the optimal compared to the other filters is that its shape and 
scale are physically and statistically well motivated such that it 
needs not be experimentally rescaled. 



In Fig. m we compare the filters studied here and in the 
literature. They are scaled in such a way as they are typically 
discussed or applied in the literature (see also the figure leg- 
end and caption for more detail). At first sight, the scales are 
surprisingly different. When the optimal filter is constructed 
including the linear matter power spectrum such as to best sup- 
press the LSS contribution, it shrinks considerably. It is reas- 

jroposed 



suring that the truncated NFW-shaped filter (THS) 
and heuristically scaled by lHennawi & Sperg el (■2005) to yield 



(36) 



best results almost exactly reproduces the optimal filter. They 
are therefore expected to perform similarly well. The optimised 
aperture-mass filter (OAPT) also peaks at fairly small angular 
scales, but shows the long tail typical for the NFW profile. The 
aperture mass has its maximum at comparatively large radii, 
explaining why the APT filter yields results most severely af- 
fected by the LSS. 
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5. Results 

5.1. Signal-to-noise maps 

We now use the above-mentioned weak-lensing estimators to 
analyse our mock catalogues of lensed galaxies. 

In practice, the integral in Eq.|2T]is replaced by a sum over 
galaxy images. Moreover, since the ellipticity e is an estimator 
for y, we can write 



AeAO) = - y £,i(0im\ei - 0\) 



(38) 



where e„(0,) is the tangential component of the observed ellip- 
ticity of the galaxy at 6i, with respect to the point 0. Similarly, 
the noise estimate in A^st is given by 



cr2(A,„)(0) = ^ 2 lf'(^«-)l'*P'(l^i 



■0\). 



(39) 



Computing „ and cr^(Aest) on a grid covering our simu- 
lated sky, we produce maps of the signal-to-noise ratio for all 
the weak lensing estimators. We use three different filter sizes 
for each estimator in order to test the stability of the results 
achieved. These have been calibrated among the different fil- 
ters to allow the optimal detection of similar objects. For the 
optimal filter, we used sizes of 1', 2' and 4'. These correspond 
to 2.75', 5.5' and 11' for the aperture mass, for the optimised 
aperture mass we used the values 5', 10' and 20' that are widely 
used in literature. 

In Fig. |5] we show examples of the signal-to-noise iso- 
contours of the weak lensing signal, superimposed on the cor- 
responding effective convergence maps of the underlying pro- 
jected matter distribution for sources at redshift Zj = 1 (left 
panels) and Zj = 2 (right panels). The iso-contours start at 
S IN - A with a step of 3. From top to bottom, the maps re- 
fer to the results obtained using the aperture mass (APT), the 
optimised aperture mass (OAPT) and the optimal filter (OPT) 
with sizes of 11', 20' and 4', respectively. The circles identify 
halos with mass M > 7 x lO''*/? 'Mq present in the field-of- 
view. The side length of each map is one degree. 

The images show that, for sources at high redshift, all three 
estimators can successfully detect the weak-lensing signal from 
clusters in the mass range considered. However, spurious de- 
tections, corresponding to high signal-to-noise peaks not asso- 
ciated with any halo, also appear. Their significance and spatial 
extent is larger in the case of the APT and th e OAPT filters. 
This confirms the results of iMaturi et al.l (l2005h . 

For lower-redshift sources, the OPT detects five out of the 
seven halos present in the field, while the APT and the OAPT 
detect substantially fewer halos. For the OPT, the number of 
spurious detections is roughly the same or slightly smaller than 
for sources at higher redshift, while it is strongly reduced for 
the APT and the OAPT. The natural explanation of these re- 
sults is that the detections with the APT and the OAPT are 
strongly contaminated by the noise from large-scale structure 
lensing, which becomes increasingly important for sources at 
higher redshift. This noise is efficiently filtered out by the OPT. 




Fig. 5. Maps of the effective convergence for sources at redshift 
= 1 (left panels) and Zs - 2 (right panels) for a region of 
simulated sky. Superimposed are the iso-contours of the signal- 
to-noise ratio of the weak-lensing signal measured with three 
estimators, namely the APT (top panels), the OAPT (middle 
panels) and the OPT (bottom panels). The iso-contours start 
from S IN - A with a step of 3. The positions of the halos 
contained in the field-of-view having mass M > 7xlO'^/!"'M0 
are identified by circles. The filter sizes are 11', 20' and 4' for 
the APT, the OAPT and the OPT, respectively. 



5.2. True and spurious detections 

In the following, we call a detection a group of pixels in the 
S IN maps above a threshold S IN ratio. Its position in the sky 
is given by the most significant pixel, i.e. that with the highest 
5 /A^ ratio. 

A true detection is obviously a detection that can be associ- 
ated with some halo in the simulation. A spurious detection is 
instead mimicked by noise, in particular by cosmic structures 
aligned along the line-of-sight. 

The association between weak-lensing detections and clus- 
ter halos is established by comparing their projected positions 
on the sky. This causes a problem, because the simulation 
boxes contain plenty of low-mass halos that are not individu- 
ally detectable through lensing but happen to be projected near 
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the line-of-sight towards a detection. Thus, spurious detections 
could easily be erroneously associated with these low-mass ha- 
los on the basis of the projected position only. 




Fig. 6. Map of the S/N ratio corresponding to a region of 3 
square degrees. The map was created using the OAPT estima- 
tor, with a filter scale of 20' and assuming a source redshift of 
Zj = 1 . The left panel shows the S /N ratio map including all 
lens planes, while the right panel shows the same map obtained 
after removing the lens plane containing the cluster responsible 
for the highest S /N peak in the left panel. 



classify the detection as true, and otherwise as spurious. We 
estimate through several checks of detections associated to the 
halos that S /N fluctuations of order 25% of the initial value are 
possible due to different properties of the noise. Thus, we set 
this limit as our threshold for discriminating between true and 
spurious detections. 

This method also shows its power when pixels identifying a 
true detection are compared with pixels associated to a spurious 
detection. This is shown in Fig. |7] The map in the left panel 
represents a true detections, while the map on the right panel 
shows a spurious detections. The maps refer to different regions 
of aS /N map created with the APT estimator with a filter size 
of ir and a source redshift of Zs = 2. As it is clearly seen, it is 
impossible a priori to distinguish which of the two is spurious. 

5.3. Statistical analysis of the detections 





Fig. 7. Maps of the S /N ratio corresponding to a region of 3 
square degrees. The maps were created with the APT estimator, 
with a filter scale of 11' and assuming a source redshift of Zs = 
2. The left panel shows a true detection, while the right panel 
shows a spurious detection. 

As pointed out earlier, we describe the lensing effect of the 
matter contained in the light cone with a stack of lens planes. 
Cluster halos are localised structures, i.e. their signal originates 
from a single lens plane. Thus, any detection should disappear 
when its plane is removed from the stack. Conversely, spurious 
detections are not caused by localised structures and should re- 
main even after removing an individual lens plane. This is il- 
lustrated by the S /N maps shown in Fig. |6l The map in the left 
panel includes all lens planes, while one plane was removed for 
the right panel. Both maps were obtained with the OAPT esti- 
mator with a filter size of 20' and a source redshift of Zs = 1 ■ 
Clearly, the highest peak in the left panel, which is in fact pro- 
duced by a massive halo, disappears in the right panel, after re- 
moving the lens plane from the stack which contains the halo. 
All other features in the left upper map remain unchanged. 

This allows us to verify the reliability of detections associ- 
ated with some halo in the catalogue. For each positive match, 
we estimate the lensing signal before and after removing the 
plane containing the candidate lensing halo from the lens-plane 
stack. If this causes a significant decrease in the S/N ratio, we 
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Fig. 8. Number of detections as a function of the S /N ratio ob- 
tained by using the APT (top panels), the OAPT (middle pan- 
els) and the OPT weak lensing estimators. Results for sources 
at redshift = 1 and Zj = 2 are shown in the left and in the 
right panels, respectively. Different line styles refer to three dif- 
ferent filter sizes. For the OPT, these are 1', 2' and 4'. They 
correspond to 2.75', 5.5', and 11' for the APT and to 5', 10' 
and 20' for the OAPT. 



In Fig. [8] we show the number of detections per square de- 
gree in S IN ratio bins, ignoring for now the distinction between 
true and spurious detections. Left and right panels refer to sim- 
ulations with sources at redshifts Zs - 1 and Zs - 2, respec- 
tively. From top to bottom, we show the results for the APT, 
the OAPT and the OPT estimators. In each panel, we use solid, 
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dashed and dotted lines to display the histograms correspond- 
ing to increasing filter sizes. 

For low source redshifts and small filter sizes, the APT 
and the OPT estimators lead to similar numbers of detections. 
Instead, for the OAPT, the number of detections is larger by 
up to a factor of two for S /N > 4. Increasing the filter size, 
the number of detections generally increases for all estimators, 
especially for large S /N ratios and in particular for the OPT. 

We notice, however, that for small S /N ratios, larger filters 
produce lower numbers of detections for the APT and for the 
OAPT. This behaviour is more evident for sources at higher 
redshifts. For example, we find that the number of detections 
with S IN - A drops by a factor of 4 for the APT and by a factor 
of ~ 7 for the OAPT, when increasing the filter size from 2.5' to 
11' and from 5' to 20', respectively. Increasing the filter size, 
the weak-lensing signal is estimated by averaging over more 
background galaxies. Thus, high S IN peaks are smoothed, and 
some detections may be suppressed. This affects mainly the de- 
tections with the APT and the OAPT filters. On the other hand, 
the OPT filter shrinks in response to the noise introduced by 
the large scale structure, largely reducing this effect compared 
to the APT and the OAPT. 

The fractions of spurious detections are shown in Fig. |9] 
Clearly, the OPT estimator performs better than the APT and 
the OAPT. For sources at redshift Zs - 1 and - 2, the fraction 
of spurious detections with the OPT is less than 20% and 30% 
at S IN ~ 4. This fraction decreases below 10% for S IN > 5 
and drops rapidly to zero for higher S IN ratios. Results are 
very stable against changes in the filter size. Conversely, the 
APT and the OAPT estimators yield similarly low fractions of 
false detections only for the smallest apertures. 

Depending on the filter shape, its size and on the source red- 
shift, aS IN threshold can be defined above which there are no 
spurious detections and thus all detections are reliable. For the 
OPT estimator, this minimal signal-to-noise ratio is between 5 
and 8. It increases above 10 for the APT and the OAPT estima- 
tors if l arge filter size s are u sed. These results agree with the re- 
su lts of Maturi et al. ( 2005h . using numerical simulations, and 
of iMaturi et al.l (120061) . regarding the analysis of the GaBoDS 
survey. 

Here, we studied the contaminations by the LSS, the intrin- 
sic ellipticity and the finite number of background galaxies all 
together. To gain an idea which of those is the main source 
for spurious detections, we used the APT with = 11' to 
analyse a catalog of galaxies with intrinsic ellipticities set to 
zero. In this case, the S/N ratio is enhanced by a factor of four 
uniformly across the whole field, but the morphology of the 
map is not affected. The same should apply to the finite num- 
ber of background sources. We thus conclude that the main 
source of spurious detections is the LSS, as already no ted by 
Reblinsky & BartelmannI (1 19991) and lWhite etaP (l2002h . 
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5.4. Sensitivity 

We shall now quantify which halo masses the weak-lensing es- 
timators are sensitive to. 



Fig. 9. Fraction of spurious detections as a function of the S IN 
ratio obtained by using the APT (top panels), the OAPT (mid- 
dle panels) and the OPT weak-lensing estimators. Results for 
sources at redshift Zs - I and Zs-^ are shown in the left and the 
right panels, respectively. Different line styles refer to three fil- 
ter sizes. For the OPT these are 1', 2' and 4'. They coiTespond 
to 2.75', 5.5', and 1 1' for the APT and to 5', 10' and 20' for the 
OAPT. 



Figure [TOl shows the lowest mass detected in each redshift 
bin. This is defined as the mean mass of the ten least massive 
halos detected in this bin. Again, results are displayed for all 
weak-lensing estimators, for different filter sizes and for two 
source redshifts. 

We note that the performance of the three filters is very 
similar for sources at redshift Zs = 1 (left panels). The OPT 
(bottom panels) is only slightly more efficient in detecting low- 
mass halos than the APT (top panels) and the OAPT (middle 
panels). The minimal mass detected depends on the lens red- 
shift. All filters allow the detection of low-mass halos more 
efficiently if these are at redshifts between 0.2 and 0.5, i.e. at 
intermediate distances between the observer and the sources. 
This obviously reflects the dependence of the geometrical lens- 
ing strength on the angular-diameter distances between the 
observer and the lens, the lens and the sources, and the ob- 
server and the sources. The lowest detected masses fall within 
~ lO^h-^Ma and ~ IO^^/i^'Mq for the OPT estimator. 

For sources at higher redshift, the region of best filter per- 
formance shifts to higher lens redshift, between 0.5 and 0.8. 
We note that due to the increasing importance of lensing by 
large-scale structures, the differences between the estimators 
are more significant. The OPT estimator allows the detection 
of halos with masses as low as < lO'^/i^^M©, almost inde- 
pendently of the filter size. Similar masses are detected with 
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Fig. 10. Minimum detected halo mass as a function of redshift 
for the APT (top panels), the OAPT (middle panels) and for the 
OPT (bottom panels) estimators. Results for sources at redshift 
Zs - ^ and Zs-2 are shown in the left and in the right panels, re- 
spectively. Different line styles refer to three filter sizes. For the 
OPT, these are 1', 2' and 4'. They correspond to 2.75', 5.5', and 
11' for the APT and to 5', 10' and 20' for the OAPT. Results 
for each redshift bin are averaged between two planes. 



the OAPT only for the smallest apertures. With the APT and 
the OAPT, the results are indeed much more sensitive to the 
filter size than with the OPT. Increasing the filter size pushes 
the detectability limit to larger masses. Again, as discussed in 
Sect. 15.31 this is due to the fact that the signal from low-mass 
halos is smeared out by averaging over an increasing number of 
galaxies entering the aperture. For example, the minimal mass 
detected with the OAPT filter at z ~ 0.8 changes by one order 
of magnitude by varying the filter scale from 5' to 20'. 

5.5. Completeness 

We now discuss the completeness of a synthetic halo catalogue 
selected by weak lensing. 

Figure [TT] shows the fraction of halos contained in the light 
cone that are detected with different weak lensing estimators 
as a function of their mass. Again, we find that the OPT fil- 
ter yields the most stable results with respect to changes in 
the filter size. This is particularly evident for sources at red- 
shift Zs - 2 (right panels), while the differences are smaller 
for Zj = 1 (left panels). As discussed earlier, the APT and the 
OAPT become less efficient in detecting low-mass halos when 
the filter size is increased. 

For the OPT estimator, the completeness reaches 100% for 
masses M > 3 x 10' and M > 2 x 10'**r'Mo for 



sources at redshift Zj = 1 and Zs - 2, respectively. For lower 
masses, the completeness drops quickly, reaching ~ 50% al- 
ready at M ~ 2 X IO'^'/i^'Mq for low-redshift sources, and at 
M ~ 7 X lO'^/z^'Mo for high-redshift sources. Similar results 
are obtained with the APT and the OAPT only for small aper- 
tures. 
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Fig. 11. Fraction of detections as a function of the halo mass. 
Each plot contains results obtained with the three filter radii 
used in this work. The panels on the left show curves for 
sources at Zs = 1, the panels on the right for sources at Zs = 2. 
From top to bottom we have the APT, the OAPT and the OPT. 



Figure [TT| gives a global view of the halos detected, regard- 
less of the their redshift. In Fig. [12] we selected three mass bins 
(M = 2.5 X lO'^Mo/Zz, M = 5 X IO'^Mq//?, M = IO'^'Mo/Zz) 
and determined the fraction of halos detected as a function of 
the redshift. 

To reduce the noise, we binned together two lens planes as 
in Fig. [To] Yet, the results are still noisy, there is much varia- 
tion for all the filters when the filter radius is changed, and the 
performance of the filters is quite similar in this respect. We 
see from the figure that the detected halos are preferentially lo- 
cated at low and moderate redshifts, due, as already said, to the 
geometry of the lensing strength. 

5.6. Comparison with the peak statistics 

The peak statistic counts peaks in convergence maps, 
e.g. obtained w ith the K aiser-Squires inversion (see 



Kaiser & Squires L 119931: iKaiser. Squires & Broadhurst 



19951) . usually smoothed with a Gaussian kernel. Even though 
they used a different set of numerical simulations, we can 
safely compare our results with the peak-statistic analysis by 
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Fig. 12. Fraction of halo detections with the APT, OAPT and 
OPT (from top to bottom) as a function of the halo redshift for 
three particular masses. The red line corresponds to a mass of 
M = 2.5 X 10^^Mo//i, the green line to M = 5 x 10^^Mo//z and 
the blue line to M = IO^'^Mq/Zz. The panels on the left show the 
results for sources at = 1 and those on the right for sources 
at = 2. From top to bottom, the filter radii are r - 5.5' (for 
APT), r = 10' (for OAPT) and r = 2' (for OPT). Results for 
each redshift bin are averaged between two planes. 



searches for dark matter concentrations in the GaBoDS survey 
dSchirmer et alll2QQ3l:lMaturi et al.Ll2QQ6b . 



To this goal, we perform a new set of ray-tracing simula- 
tions, where a realistic redshift distribution of the sources is as- 
sumed. In particular, we draw the sources from the probability 
distribution function 



P(z)=iVexp[-(z/zo)^], 
where is chosen such that 



Jo 



P(z)dz = 1 



(40) 



(41) 



We adapt P(z) to the redshift distribution of the sources in 
th e GaBoDS survey by setting zo = 0.4 and /3 = 1.5 



(ISchirmer et all 120031) . In order to mimic the number den- 
sity of galaxies in the GaBoDS observations, we assume rig = 
lOarcmin"'. 

By repeating the same analysis outlin ed above, we find re - 
sults that are compatible with the results of Maturi et al. ( 2006h . 
In particular, the number of detections with S /N = 3.5 per 
square degree in our GaBoDS simulations (in GaBoDS data) 
ai-e ^ 5 (^ 4) for the OPT with r = 2', ^ 3 (^ 3) for the OAPT 
with r = 10' and ^ 1.5 2) for the APT with r = 5.5' (r = 4') 
respectively. A comparison between the detections with differ- 
ent weak lensing estimators is shown in Fig. [T3] 
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Hamana et al.l (120041) . whose Gaussian kernel has a FWHM of 
1 arcmin. 

Fix ing a detection threshold of S /N > 4 (5), 'Hamana et al.l 
(120041) found N ~ 6 (2.5) detections per square degree, 60% 
(76%) of which correspond to real haloes with masses larger 
thanlO'^/i^'Mo.In our simulati ons, with th e same S/N thresh- 
old and the optimal filter by Maturi et al.l ( l2005b . we found 

^ 10 (7), with an efficiency in detecting real haloes of 
85% (95%). For halos with masses M > 2 x IO'^/z^'Mq 
(M ^ 1x10'^ /i 'Mg). the Hamana et al] (|2004i) sample is com- 
plete at the 70% (50%) level, which is virtually identical to the 
completeness of 70% (50 - 60%) achieved with the optimal 
filter. 
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6. Comparison with observations 

The results outlined above show interesting differences be- 
tween the performances of the filter functions. The discrepan- 
cies are particularly significant for high-redshift sources, in- 
dicating that the noise due to the LSS should become impor- 
tant only for deep observations. We can now attempt a quick 
comparison of our simulations with the observational results 
existing in the literature. In particular, we focus here on the 



Fig. 13. Total number of detections per square degree (left 
panels) and fraction of spurious detections (right panels) 
for sources distributed in redshift as in the GaBoDS survey 
(ISchirmer et al.Ll2003h . From top to bottom, we show the APT 
(for r = 2.75', r = 5.5' and r = 1 1'), the OAPT (r = 5', r = 10' 
and r = 20') and the OPT (r = 1', r = 2' and r = 4'). 

The fraction of spurious detections is large for all filters, 
but it is generally smaller for the OAPT and the OPT. As ex- 
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pected, the OAPT and the OPT estimators have similar perfor- 
mances, because of the small density of background galaxies. 
Indeed, the noise due to the intrinsic shape of the sources is 
dominant with respect to that due to the LSS and thus, accord- 
ing to Equation ([34]), the two filter functions have a very similar 
shape. 

A more detailed statistical analysis, including simulations 
representative of different current surveys (both ground- and 
space-based), will be presented in a future paper. 

7. Conclusions 

We studied the performance of dark-matter halo detection with 
three different linear filters for their weak-lensing signal, the 
aperture mass (APT), the optimised aperture mass (OAPT), and 
a filter optimised for distinguishing halo signals from spurious 
signals caused by the large-scale structure (OPT). In particular, 
we addressed the questions how the halo selection function de- 
pends on mass and redshift, how the number of detected halos 
and of spurious detections depends on parameters of the obser- 
vation, and how the filters compare. 

To this end, we used a large A^-body simulation, identified 
the halos in it and used multiple lens-plane theory to deter- 
mine the lensing properties along a fine grid of light rays traced 
within a cone from the observer to the source redshift. Halos 
were then detected as peaks in the filtered cosmic-shear maps. 
By comparison with the known halo catalog, spurious peaks 
could be distinguished from those caused by real halos. 

Our main results are as follows: 

- We confirm that a mong those te s ted, th e optimised filter 
(OPT) proposed by Maturi et al. ( 2005h performs best in 
the sense that its results are least sensitive to changes in the 
angular filter scale, it produces the least number of spurious 
detections, and it has the lowest mass limit for halo detec- 
tion (cf. Figs.|9l[l0]and[n]). 

- With the OPT filter, the fraction of spurious detections is 
typically < 1 0% for a signal-to-noise threshold of S /N ~ 5. 
It increases with source redshift due to the larger contami- 
nation by large-scale structure lensing (cf. Fig.|9]l. 

- The number of halos detected per square degree by the OPT 
filter with 5 /A^ > 5 is a few if the sources are at redshift 

= 1, and ~ 20 for - 2 (cf. Fig.|8]l. 

- The minimum detectable halo mass starts at a few times 
lO'^/i^'Mo at redshifts ~ 0.1, drops to ~ 10'^/;"' Mq 
near the optimal lensing redshift and increases towards 
~ lO'"^/! ' Mq approaching the source redshift (cf. Fig.fTOll. 

- The fraction of halos detected reaches ~ 50% at ~ 2 x 
lO^'^h-^ Mq and 100% at ~ 4.5 x lO'"*/!"' Mq with sources 
at Zs - 1 . With more distant sources at Zs - 2, half of the 
halos with ~ 7 x lO^^h'^ Mq are found, and all halos above 
~ 3 X 10'4/i-i Mo (cf. Fig.in)- 

- Adapting par a meters to the GaBoDS survey 



jSchirmer et al.L 120031) . and distributing sources in 
redshift, our simulation yields a number of significant 
detections per square degree which is in good agreement 
with what was found applying the OPT filter to the real 
GaBoDS data dMaturiet al.-.2006,) . 



Thus, the OPT filter, optimised for suppressing contam- 
inations by large-scale structures, allows the reliable detec- 
tion of dark-matter halos with masses exceeding a few times 
lO'-' /i"' Mq with a low contamination by spurious detections. 
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